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Abstract A variational solution to the rovibrational problem in curvilinear vibrational co- 
ordinates has been implemented and used to investigate the nuclear motions in several linear 
triatomic molecules, like HCN, OCS, and HCP. The dependence of the rovibrational energy 
levels on the rotational quantum numbers and the /-doubling has been studied. Two approx- 
imations to the rovibrational Hamiltonian have been examined, depending on the level of 
truncation of the potential energy operator. It turns out that truncation after the fifth order 
in the potential is sufficient to produce vibrational energies of high accuracy. An interest- 
ing feature of the present formulation of the problem in terms of the curvilinear vibrational 
coordinates is the explanation for the /-doubling of the rovibrational levels, which in this 
picture is interpreted as the result of the inequivalency of the average rotational constants 
in mutually perpendicular planes, rather than as the effect of the Coriolis-type interactions 
between the vibrational and rotational motions. The present theoretical results are compared 
with the available experimental data from high-resolution spectroscopy, as well as with other 
ab initio calculations. 

Keywords rovibrational spectra • curvilinear vibrational coordinates • anharmonicity • 
/-doubling • HCN • OCS • HCP 



1 Introduction 

A proper description of the interactions between the rotational and vibrational motions is an 
important issue in molecular spectroscopy and the first attempts to find a theoretical solu- 
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tion date back to the beginning of the previous century. Even a short overview of the most 
important methods is beyond the capacities of a single research paper, so we therefore refer 
the reader to a recent tutorial overview of the existing approaches [1], and restrict ourselves 
to a necessarily incomplete review of the literature referring to successful implementations 
of the methods directly related to the topic of our paper. 

The choice of the variational method has several advantages over the perturbational the- 
ory alternative. Although in many cases perturbation theory allows for a physically sound 
interpretation of the results, as is the case with the rovibrational energy levels, its possi- 
ble drawbacks involve among other things the convergence problems of the perturbation 
series, which do not appear if the variational alternative is used. The variational method is 
therefore much more robust, especially if many states are to be treated at once and if Fermi 
resonances are expected. One of the earliest attempts to solve the vibrational problem with 
the variational method was reported in Ref. [2] utilizing the Watson Hamiltonian [3]. Still 
other applications of this type can be found, e.g., in Refs. [4, 5], A solution to the full rovi- 
brational problem for a linear molecule based on the isomorphic Watson Hamiltonian in the 
linear normal coordinates with a basis set spanned by the 2D harmonic oscillator eigenfunc- 
tions in a polar form is in being sought in our laboratory [6]. In Ref. [7] a method to solve 
the variational rovibrational problem for linear triatomic molecules has been reported, which 
was then developed in Refs. [8, 9]. However, the methods mentioned above employ internal 
vibrational coordinates with the exact vibrational Hamiltonian, which make them difficult to 
extend to systems with more than three atoms. The mixed variational-perturbative approach 
used by Suzuki [10-12] should also be mentioned in this context. 

The well-known variational approach to triatomic molecules was implemented by Ten- 
nyson etal. [13, 14] within the DVR3D software package, which allows the solving of the full 
rovibrational problem. An example application of DVR3D for the HCN and HNC molecules 
can be found in Ref. [15]. However, in the method of Tennyson et al. the Jacobi or Radau 
coordinates are used for the exact (within the Bohr-Oppenheimer approximation) Hamilto- 
nian, which makes it difficult to extend to molecules with a larger number of atoms. It is 
also worth mentioning the Hougen-Bunker- Johns approach for nonrigid molecules [16, 17] 
implemented in the TROVE package [18]. The illustrative application of this approach can be 
found in Ref. [19] for HDO and Ref. [20] for HSOH. However, in the Hougen-Bunker- Johns 
approach linear molecules should be treated in a special way [21], and we are not aware of 
any implementation of this method for this particular (linear) case. 

The method presented in this paper is more universal and can be applied to larger 
molecules, both linear and nonlinear. It was proposed for the first time by Pavlyuchko [22] 
and is in principle limited only by the computational costs of the calculation of the potential 
energy surface (PES) and by the dimensions of the largest block of the resulting rovibra- 
tional Hamiltonian matrix. Some attempts to implement this method for linear molecules 
have been undertaken recently [23], but up to now no results obtained with this approach are 
available in the literature. It is worth noting that many features of the rovibrational problem 
depend on the way in which the vibrational and rotational coordinates are introduced. The 
interpretation of the /-doubling is one of the prominent examples [24, 25]. 

The molecules selected for the first test of the new implementation of the rovibra- 
tional motion problem are three triatomic molecules: HCN, HCP, and OCS. The first two 
molecules contain a light hydrogen atom and show a large amplitude vibration, while OCS 
is a rigid molecule with a low amplitude of the bending mode. Additionally, for the HCN 
and HCP pair, the effects of introducing a heavier atom belonging to the same group can be 
examined. 
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All three selected molecules are of considerable astrophysical and astrochemical inter- 
est. Among them HCN is undoubtefully the most widespread one, and numerous experi- 
mental and theoretical studies have been undertaken regarding this molecule (see e.g., Refs. 
[26-30]). It is also one of the most prominent interstellar and circumstellar molecule. The 
degenerate bending mode is an important feature of HCN because of the large amplitude of 
motion of the light hydrogen atom and also because the bending motion is a direct pathway 
to the HNC isomer. 

Theoretical and experimental data for the OCS molecule can be found in Refs. [26, 31— 
33]. Spectroscopic interest in the OCS molecule stems from the fact that this molecule is 
present in the Venus atmosphere along with the more abundant carbon dioxide [34]. Exper- 
imental high-resolution rovibrational spectra have been recorded and assigned e.g., in Refs. 
[35^10], in some cases also for less common isotopomers of OCS. One of the first theoreti- 
cal analyses of the OCS spectra utilized the semiclassical approach based on the integration 
of selected trajectories [41]. Martin et al. [42] used the quartic force field obtained from 
the coupled cluster calculations including single, double, and noniterative triple excitations, 
CCSD(T) [43], with the cc-pVQZ basis set to find the vibrational levels and some rotational 
constants for the OCS, CS, H2S and CS2 molecules by the second-order perturbation theory. 
They found for the carbonyl sulfide molecule a generally good agreement with the available 
experimental data [44], although some Fermi resonances had to be accounted for in order to 
get nearer to the experimental energy levels. For instance, the C=0 stretching fundamental 
band (Vi) interacts with the 02° 1 combination band. Peterson et al. [45] obtained a stretch- 
ing PES for the carbonyl sulfide molecule with the fourth-order M0ller-Plesset perturbation 
theory (MP4) and the size-consistency corrected configuration interaction method limited 
to single and double excitations (CISD) and used them to calculate the stretching bands of 
OCS as well as the rotational and /-doubling constants. Both stretching and bending modes 
have been examined by Pak and Woods [46], who used the CCSD(T) method to obtain the 
PES of OCS and the second-order perturbation theory for the spectroscopic constants. More 
recently, Xie et al. [47] devised a PES capable of reproducing high energy vibrations by 
refining the force field constants of Pak and Woods [46] through the fitting to experimental 
vibrational levels up to 8000 cm~ 1 . 

The methylidynephosphine molecule HCP is an unstable species which decomposes af- 
ter several hours in room temperature conditions [48], but it belongs to a small number of 
phosporus-containing molecules present in the interstellar media [49] and is therefore of 
considerable astrochemical interest. From the chemist's point of view this molecule is a rare 
example of the existence of the unusual C=P triple bond. Interest in the spectroscopic char- 
acteristics of HCP grew in the early nineties of the last century due to speculations about the 
presence of this molecule in some planetary atmospheres [50]. For its detection in such en- 
vironments a detailed knowledge of the laboratory rotation-vibration spectrum is required. 
Experimental examination of the HCP spectra has been reported e.g. in Refs. [51-55]. The 
quality of these measurements ranges from the earliest low-resolution spectra to the most 
recent high-resolution ones. The theoretical work on this molecule usually included the cal- 
culation of the PES, which in some cases was then utilized to reproduce vibrational levels 
and rotational constants [52, 54, 56, 57]. Koput and Carter [56] used the variational method 
for the vibrational problem followed by the perturbation theory to describe the coupling be- 
tween rotations and vibrations. Specifically, perturbation theory was applied to calculate the 
rotational constants with the PES calculated at the CCSD(T) and the second-order M0ller- 
Plesset perturbation theory (MP2) levels of theory. Koput and Carter [56] concluded that 
the MP2 level is not sufficient to reproduce the rovibrational spectra, while the CCSD(T) 
method is capable of doing so provided that large enough orbital basis is used in the calcu- 
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lations. The basis set effect is particularly important for the vibrational energy levels, while 
the anharmonic constants are reproduced quite well in the smallest cc-pVTZ basis. Koput 
and Carter [56] also tested the influence of the correlation contributions of the core elec- 
trons to the vibrational and rotational parameters. They concluded that the core electrons' 
correlation shifts the vibrational levels to higher energies and has little influence on other 
constants. This line of study has been continued by Beck et al. [57] who calculated the po- 
tential energy surface with the multireference configuration interaction method limited to 
single and double excitations (MRCI) and used it to calculate high vibrational states and the 
first rotational constant B. Among other theoretical works on HCP one should mention the 
study of polyads of highly excited vibrational states with the Fermi resonance Hamiltonian 
[58]. 



In this paper the implementation of the Rayleigh-Ritz variational solution for the rovi- 
brational problem with the Hamiltonian expressed in the curvilinear vibrational coordinates 
and the Euler angles is reported and applied to a number of small triatomic molecules. The 
simplicity of the proposed method follows from the polynomial representation of the po- 
tential energy surface and the one-dimentional (ID) harmonic oscillator functions in the 
basis set. In contrast to most of the other approaches employing the curvilinear coordinates, 
here the variational method is used for the full rovibrational problem, and not just for its 
vibrational part, which allows for the larger flexibility of the method compared to the pertur- 
bation theory. With the present implementation we calculated selected rovibrational levels 
of HCN, HCP, and OCS up to 13000 cm"'. For higher energies the polynomial represen- 
tation of the potential energy operator becomes inappropriate and one needs a Morse-like 
potential and Morse oscillator functions in the basis set. In this case the matrix elements of 
the Hamiltonian become much more complicated [59, 60]. However, we believe that the cur- 
rent implementation based on the polynomial expansion is sufficient to obtain and interpret 
many low-energy rovibrational spectra. In our investigation we paid a special attention to 
the rovibrational interactions, caused by the dependence of the inverse inertia matrix on the 
vibrational coordinates. The rovibrational levels were obtained up to and including / = 15, 
then fitted to the known functional form including rotational centrifugal constants and con- 
stants of the /-resonance and, finally, compared with the available experimental data. 



The plan of this paper is as follows. In section 2 we present the mathematical approach 
to solve the problem of nuclear motion in polyatomic molecules based on the Hamiltonian 
written in the curvilinear vibrational coordinates and show some advantages of this Hamil- 
tonian over a Hamiltonian expressed in terms of linear vibrational coordinates. In section 3 
we report the details of the Rayleigh-Ritz variational procedure, discuss the choice of the 
basis set functions, and the way the potential energy function should be treated. In section 
4 the analytical fitting of the rovibrational energy levels, which are widely used to represent 
the experimental data, are discussed. The same analytical expressions are utilized to fit the 
computed energy levels. Section 5 gives an explanation of the /-doubling effect when the 
curvilinear coordinates are used. In section 6 we give some computational details of the po- 
tential energy calculations and of the matrix diagonalization in the variational procedure. In 
section 7 we compare our results with previous theoretical calculations and with the data 
derived from high-resolution spectroscopic experiments. Finally, section 8 concludes our 
paper. 



5 



2 Molecular Hamiltonian 

Several ways of dealing with the molecular Hamiltonian describing rovibrational motions 
have been described in the literature, depending on whether linear or curvilinear vibrational 
coordinates are used in the vibrational part. In this section the theory in curvilinear coordi- 
nates will be presented. We start the theory section by introducing the curvilinear vibrational 
coordinates Q s [59,61] defined as linear combinations of the natural (internal) vibrational 
coordinates q, : 

The matrix L describes the mathematical form of the s-mode vibration [62, 63]. The matrices 
L and L p are related to one another by the following expression: 

(L") r L = I, (2) 

where I stands for the unit matrix. 

The coordinates Q s are also known as curvilinear vibrational coordinates, since they 
are related nonlinearly to the Cartesian coordinates of the atoms. For linear molecules the 
natural vibrational coordinates qi are of two types, corresponding to the bond stretching and 
to the linear angle bending modes. The latter mode, which is doubly degenerate, is described 
by two coordinates in the mutually perpendicular planes [63]. The kinetic energy operator 
for a polyatomic molecule can be written as [64, 65]: 

f#lji** gij{q)g -U2d y 4 

2 dqi ( )cjj 

The vibrational part of the Hamiltonian can be reduced by performing the differentiation 
[59, 65]: 

where gjj (q) denote elements of the kinetic matrix and /3 (q) is a nondifferential kinematic 
operator [59], also referred to in the literature to the pseudo-potential [66]. The components 
of the g-tensor are well known and can found in Ref. [64]. Some features of g for linear 
molecules are explained in Refs. [67, 68]. 

The exact dependence of the g-matrix on the vibrational coordinates is rather compli- 
cated, therefore some approximations should be employed in Eq. (4) in order to obtain 
tractable formulas. Different approximate treatments of H v in the curvilinear coordinates 
have been proposed in the literature [59]. In this work we chose to expand the elements 
gij(q) in the Taylor series: 

«/(*) = «,(o)+l(%^) ft+Lf^) **+•••• <» 

k\ d 1kJo ki\ d 1k<H Jo 

If additionally the pseudopotential j3 (q) is neglected, the following expression for the vibra- 
tional Hamiltonian is obtained: 

Ti* d 2 h 2 ^ d d fc* d d 
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where the kinematic coefficients Xspr and Xsprt in Eq. (6) are defined by the formulas:, 



-■s(^).*^ 



ijkl 

It should be noted that for linear molecules the matrix of the L^p coefficients is composed 
of two submatrices, corresponding to stretching and degenerate bending modes with no 
cross-terms. This explains why the only non-vanishing components of % are Xsps and Xspp's' 
the letters s and p refering to the bending and stretching coordinates, respectively. 

The total Hamiltonian describing the rovibrational motions of a polyatomic molecule 
can be expressed through the curvilinear vibrational coordinates and the Euler angles (0 , <j> , £ ) 
which describe the rotation of the axes of the equilibrium inertia tensor with respect to the 
laboratory frame [22]: 

h 2 

H = H v + -^^i a p(Q)J a Jp, (8) 

where P- a p{Q) denote elements of the inverse inertia tensor matrix. To obtain a working 
expression for Eq. (8), we expand pi (Q) at the equilibrium point in the curvilinear vibrational 
coordinates through the second order: 

a* (fi) = ^ (o) + E & + - 2 1 ( J &&+■••■ (9) 

In order to find the derivatives appearing in Eq. (9) the expansion coefficients of the inertia 
tensor in the linear normal coordinates g' r are used: 

U =%p+Zal a P ) Qi+Z4f ) QiQi+.... (10) 

r rs 

Explicit expressions for the coefficients tf* = (^) q and J&> = are well 

expressed through the al"^ and terms, cf. Ref. [59] for the details. 

We also have to note that linear molecules have 3N — 5 vibrational degrees of freedom 
and the inertia tensor l a n has only two main axes, x and y. The vibrational modes can be clas- 
sified as TV — 2 doubly-degenerate bending modes and N — 1 stretching modes. Degenerate 
modes are described by two coordinates in the perpendicular planes [63]. The internal bend- 
ing coordinate is introduced as sin0 fl , where 6 a is the angle of the distortion from linearity. 
In our ab initio calculations we used 9 a instead of sin a , because one can approximately set 
sinjc w x for small values of x. 

The rovibrational Hamiltonian expressed through the curvilinear vibrational coordi- 
nates, Eq. (4), does not contain Coriolis-type rovibrational interaction terms, i.e. terms of 
the type J a pp and p a Jp, which are present in the rovibrational Watson Hamiltonian derived 
in Refs. [3,70]. The Hamiltonian given by Eq. (8) contains just the centrifugal interaction 
terms since p: a p (Q) depends on the vibrational coordinates only. The other advantage of the 
latter Hamiltonian form is that there is no need for the special treatment of linear molecules 
as in the case of the Watson Hamiltonian. 





known in the literature [69]. The coefficients ( "^q^ 1 I and I " ^QdQ. ' I can * nen ^ e 
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3 The variational solution to the rovibrational problem 

In order to find the rovibrational energy levels we have to solve the rovibrational problem 
described by the Hamiltonian given in Eq. (8), i.e. to solve the following Schrodinger equa- 
tion: 

^(e,M,S) = #*9i(G,M,§), (ii) 

where k denotes a super-index containing quantum numbers corresponding to the vibra- 
tional (n) and the rotational (J,M,K,p) degrees of freedom, i.e. k = {n,J,M,K,p}. To solve 
the eigenproblem given by Eq. (11) we follow the algorithm described in Ref. [22], where 
the Rayleigh-Ritz variational method is applied. In this algorithm the wave function is 
expanded in a basis of known functions (pi(Q, 0, 8, |): 

«ftG3,M,§) = LQm(<2,0, 0,1), (12) 

where the basis function <Pi(Q,$,6,^) is a product of a vibrational function yf n (Q) and 
rotational function a J MK (§, & ,£,): 

<p l (Q,<l>,0,Z) = y n {Q)ai jK (<l>,0,Z). (13) 

The functions a J MK (§, 6,%) are solutions to the Schrodinger equation for the symmetric top 

a J MK {^e,^) = ^^D^,e^), (14) 

where Dy K ($ ,0 ,£;) are the Wigner D-functions [71]. The dimension of the basis set {<jP/} 
can be minimized if the vibrational functions y n {Q) are variational solutions to the multidi- 
mensional vibrational Schrodinger equation: 

H v Wn(Q)=E„Wn(Q)- (15) 

In Eq. (15) n stands for several vibrational indices, e.g. for a linear triatomic molecule 
n = {rti,n2a,n2b 1 n^}. Here nia^n^b are the indices of the basis set functions describing the 
motion in the perpendicular planes of the degenerate bending mode. 

The vibrational problem, Eq. (15), is solved in a basis set constructed from the products 
of the ID harmonic oscillator eigenfunctions (j>i(Q) for each vibrational degree of freedom: 

= Xn s H ns {Qs) exp[-g?/2] , 0, e] - oo, oo[, (16) 

where H„(Q S ) is the Hermite polynomial of the nth order and N„ = (2"n!y / 7r) -1 / 2 is the nor- 

,.1/4 

malization constant. We have also introduced Q r , a dimensionless coordinate, Q r = Q r ^j%, 
k r is the harmonic force constant for the rth vibrational mode, and h the Planck constant in 
the appropriate unit system. In the particular case of a linear triatomic molecule the multidi- 
mensional vibrational wave function is expanded as: 

0«i0«2«0"26 , fVv (I 7 ) 

n 

The last remaining issue for the pure vibrational problem is the representation of the 
potential energy operator V. Usually the Taylor expansion around the equilibrium point 
truncated after the lowest few terms is used in practical calculations. The Taylor expansions 
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of the potential energy function in both the internal and curvilinear coordinates are listed 
below: 

= ^Ylfu^^j + ^^fu^^j c Jk + ^T,fijki<li c lj<lkqi + ---, (18) 

ij ijk ijkl 

V(Q) = \Y* F *pQsQp + \L F s P rQiQjQk + ^7 E FspnQsQpQrQt + ■■■■ (19) 
sp spr sprt 

With the help of Eqs. (18) and (19) the expansion coefficients /, (for the internal coordi- 
nates) and Fj,„ (for the curvilinear vibrational coordinates) are implicitly defined. The F;... 
coefficients are obtained from the fi„ coefficients, which in turn are calculated by an ab 
initio method. The corresponding transformation formulas for the curvilinear vibrational 
coordinates are given by the following equations [63]: 

Fip = ^f'jLisLjp, 

ij 

Fspr fijkLjsLjpLjf , 

ijk 

Fsprl — ^ fijklLi s L jpL j r L j t . (20) 

ijkl 

This transformation is simpler than the analogous expression for the linear vibrational coor- 
dinates [72] since no higher-order components of the L tensor appear in Eq. (20). Of course, 
in the harmonic approximation both approaches (linear and curvilinear) are fully equivalent. 

In a finite basis of the functions <p;, the solution of the rovibrational problem, Eq. (11), is 
reduced to the diagonalization of the Hermitian matrix with the following matrix elements: 

H nK ,n>K> = {Vn{Q)a J MjK {<t>,e,t; ) \H\ w {Q)a J MjK , (<j), B,$ )). (21) 

This general formula can greatly be simplified if (/) the structure of the rovibrational Hamil- 
tonian, Eq. (8), (if) the properties of the vibrational basis functions, and (Hi) the known 
selection rules for the Wigner D-functions are all taken into account. After some algebra 
[22], the following equation for the element of the Hamiltonian matrix is obtained: 

H„K.n'K' = E n 8, m > 8 KK i + — £ \l™ ^ (J,K\J a Jp\J,K') , (22) 
1 ap 

where the p.'™p quantities are defined as: 

p.$ = (n\^ ap (Q)\n') = J Wn(Q)^(QW,AQ)dQ. (23) 

It should be stressed that p.'™p are the coefficients that determine the rovibrational interac- 
tions. The diagonal terms in Eq. (23) can be interpreted as effective rotational constants for 
the nth vibrational state, while the off-diagonal terms are responsible for the interactions of 
different vibrational energy levels with rotational motion. All coefficients appearing in the 
equations above can be found in Ref. [22]. 

For 7 = the diagonalization of the Hamiltonian matrix gives pure vibrational energy 
levels. For states with / > rovibrational interactions appear. Their physical nature is cen- 
trifugal since the elements of the matrix fi(Q) from Eq. (8) depend on the vibrational coor- 
dinates Q. Obviously, neglecting of the dependence of the matrix fi(Q) on the vibrational 
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coordinates by taking only the first term in the expansion given by Eq. (9), makes the ma- 
trix p., Eq. (23), diagonal. In this case the rovibrational problem can be separated into pure 
vibrational and pure rotational problems. 

4 Constants of rovibrational interactions for linear triatomic molecules 

The formulas for the constants of rovibrational interactions and for /-doubling are well 
known in the literature. However, for the sake of completeness and to introduce the no- 
tation we repeat them here. The rovibrational levels for the linear top have been derived in 
Refs. [73, 74] and are taken in the present form from Refs. [28-30, 75]. 

The rotational energies for a given vibration, as obtained from Eq. (11), can be expressed 
in the following form: 

E r =E? + E?- d , (24) 

where Ef and Z?"~ d denote diagonal and non-diagonal parts of the expression, respectively. 
The diagonal part is given by the following formula: 

E? = B v [J (J + 1 ) - K 2 } - D v [J {J + 1 ) - K 2 } 2 + H v [J {J + 1 ) - K 2 } 3 , (25) 

where v denotes the vibrational quantum numbers, v = {Vj , v¥ , V3}, B v is the rotational 
constant, and D v and H v are the quadratic and cubic constants of the centrifugal distortion, 
respectively. The non-diagonal term was fitted to the following expression [30]: 

(Vi,V 2 ,V3,/,7|^|Vi,V 2 ,V 3 ,/±2,7) = ±^[q v -q VJ J(J+l)+q v jjJ 2 (J+l) 2 + q,(l±\) 2 } 

x [(v 2 ±/)(v 2 ±Z + 2)(7(7+l)-/(/±l))(7(7+l)-(/±l)(/±2))] 1/2 . (26) 

Here, ,3^ stands for some effective perturbation Hamiltonian used in Ref. [30]. In the par- 
ticular case of the (01^0) vibrational states we have: 

E?- d = ±i[ 9v -q VJ J(J+ l)+q VJ jJ 2 (J+ 1) 2 }J(J+ 1). (27) 

More formulas for the vibrational states (02°0), (02 ±2 0) and (Ov^'O) can be derived from 
Eq. (26). The nonzero constants D v ,H v ,q v , q v j, and q v jj result from the dependence of fi a p 
on Q, cf. Eq. (9). It should be noted that the variational approach allows the calculating of 
the dependence of D and H on J, while perturbation theory allows us to find these constants 
only for the ground vibrational state [2, 56]. 

5 Symmetry considerations and /-doubling 

Triatomic linear molecules have three normal modes: two stretching vibrations usually de- 
noted by Vi and V3 and one degenerate bending vibration denoted by v 2 . Depending on 
the presence or absence of a center of symmetry, linear molecules have either or C xv 
symmetry. For noncentrosymmetric molecules two stretching vibrations belong to the E + 
representation, while the bending has J7 symmetry [16]. The stretching vibrations in this 
case can usually be approximately attributed to a particular bond, e.g. for the HCN molecule 
Vi corresponds approximately to the C-H stretch and V3 to the C=N stretch. 
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A few words of explanation are due as far as the nature of the /-doubling is concerned. 
Commonly [3, 24], when the rovibrational Watson Hamiltonian for linear molecules is used 
[76], this effect is described by the Coriolis interactions of the rotational and vibrational 
motions of the molecule. The explanation of the /-doubling when the Hamiltonian is ex- 
pressed in curvilinear vibrational coordinates is different since, as has already been pointed 
out in Section 2, there is no Coriolis coupling in Eq. (8). In the current approach this effect 
is explained by the fact that for degenerate molecular deformations in two directions the 
average values of the inertia tensor are different. As a consequence, for a given vibrational 
state rotational states with the same absolute value, but with a different orientation of the 
projection of angular momentum on the symmetry axis K, have different values of the rota- 
tional constants p.'™p and therefore also different rotational energy levels. These two possible 
explanations for the /-doubling have been pointed out by Herzberg [24]. 



6 Computational details 

Force constants corresponding to the potential energy function V(q), Eq. (18), were calcu- 
lated at the CCSD(T) level with the MOLPRO suite of programs [77]. The cc-pVQZ basis 
set was used for all atoms. Only valence electrons were correlated in the CCSD(T) calcu- 
lations. Since no analytical gradients are available in MOLPRO for the CCSD(T) method, 
finite-difference formulas were applied in order to obtain the appropriate force constants. 
This technique naturally limits the maximum order of the force constants that can be calcu- 
lated with sufficient precision. In order to check the accuracy of the numerical procedure, 
both 5- and 7-point formulas were used for one-variable derivatives and multi-point formu- 
las for mixed derivatives. Then, only the converged digits were utilized in the subsequent 
calculations. 

The formulas for some exemplary finite-difference fourth derivatives used in this work 
are given below, where it is assumed that the potential energy function in the internal coor- 
dinates is given by V(n, >~2, Q a )- In the derivative calculations the distance increments h,hi, 
and ha were set equal to 0.01 A and the angle increment a to 1°. 

K^\ = (6V(0,0,0)-4(V(/i,0,0)+y(-/2,0,0)) + (V(2/i,0,0)+y(-2/2,0,0)))//i 4 
+ 0{h 6 ), (28) 



d 4 V\ 



dr\ to 



56V(0,0,0) - 39(y(/!,0,0) +V{-h,0,0)) 
+ 12(7(2/1,0,0) +y (-2ft, 0,0)) - (V (3/1,0,0) + V{-3h,0,0))"j/{6h 4 ) + ^(/i 8 ), (29) 



d 4 V 



I = (v(h h hz,a)-V(hi,-h2,d)+V(h u -h2,-a)+V(-h h h2,a) 
drfdr2od a 10 V 

-V(-h u -h2,a)-V(-h u h 2 ,-a)+V(-h u -h 2 ,-a)+2(V{0,h 2 ,a)-V{0,-h2,a) 
- V{0,h 2 ,-a) + y(0, -h 2 , -a)]\/{Ah}h 2 a) + ^(hjhla 2 ), (30) 



1 1 



d 4 V 



(l6(V{hM,a)-V(hi,-h 2 ,a)+V(hi,-h2,-a)+V(-h,h 2 ,a) 



drfdr 2 d6 a 

-V(-h h -h2,a)-V(-huh2,-a)+V(-hu-h 2 ,-a)) + 30(-V(0,h 2 ,a) 
+ V{0,-h 2 ,a)+V(0,h 2j -a)-V(0,-h 2 ,-a)) -V(2h h h 2 ,a) +V(2h h -h 2 ,a) 
+ V(2hi ,h 2 ,a)- V(2hi , -h 2 . -a) - V(-2h\ , h 2 , -a) + V{-2d x , -h 2 ,a) 

+ V(-2hi,h 2 ,-a)-V(-2h h -h 2 ,-a))/(48hjh 2 a) + 0(hjhja 2 ). (31) 



In order to estimate the accuracy of the Taylor expansion for V, the series truncated af- 
ter the fourth and the fifth orders have been used in the production calculations. They are 
referred to as Approximation 1 and Approximation 2, respectively. Application of Approxi- 
mation 3, i.e. the series truncated after the sixth order, is unfortunately not possible because 
of the limited accuracy of the numerical differentiation. Since the calculations reveal that for 
molecules involving hydrogen atoms the contributions from the fifth order are significant, 
and hence the sixth-order terms should possibly be included, we expect somewhat larger 
errors in the rovibrational energy levels for the HCN and HCP molecules. 

The rovibrational energies were calculated with a program written by one of our team 
(L.S.). The masses of the most common isotopes were used in the calculations of the rovi- 
brational spectrum. In all approximations we used the basis set containing 10000 harmonic 
functions, i.e. 10 basis set functions for each vibrational degree of freedom, which was 
enough for the convergence of at least the 200 lowest eigenvalues. The Hamiltonian ma- 
trix consists of the sub-blocks pertaining to each K = —J, J. However, these subblocks are 
not independent, since they are connected through the dependence of fi a p on the coordi- 
nates Q, and the non-diagonal terms become responsible for the rovibrational interactions. 
The Davidson algorithm [78] was used to find the lowest eigenvalues and the corresponding 
eigenvectors of the rovibrational Hamiltonian. The rovibrational energy levels were then as- 
signed by examining the largest weight of the harmonic part for the vibrational levels and 
the largest weight of the rigid-rotator eigenfunction. In this manual procedure we were able 
to localize the first 30 vibrational levels. An unambiguous identification for higher energy 
values turned out to be impossible because of the existence of anharmonic eigenfunctions 
almost equally distributed between some harmonic eigenfunctions. Special care has to be 
taken to ensure proper description of the vibration motion within the C-H bond. Because of 
the presence of the light hydrogen atom this mode is floppy and highly anharmonic. In order 
to improve this description the "extended" harmonic basis functions have been added to the 
basis for for the Vi -mode, corresponding to the C-H stretch, utilizing the scaled exponents 
ak s instead of the harmonic force constants k s in Eq. (16), where a was chosen from the 
interval (0.5-1.0) to enable better convergence. 



7 Results and discussion 

In Table 1 equilibrium distances calculated in this work for the molecules under study are 
presented and compared with the experimental values. An inspection of this Table shows 
that for all molecules the difference between the experimental and the ab initio distances 
between the carbon atom with the hydrogen atom is small, about 0.002 A. Somewhat larger 
differences, 0.007 and 0.008 A, are observed for the second distance if the third-row atom is 
connected to the carbon atom. Most probably, this is a result of neglecting the correlation of 
the core electrons. However, as pointed out in Ref. [42], such a small discrepancy has only 
a minor effect on the calculated rovibrational constants. Therefore we decided to proceed 
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with the calculations of the CCSD(T) force constants, Eq. (18), with only valence electrons 
correlated. 

For the HCN molecule the importance of various terms in Eqs. (5)-(6) on the com- 
puted vibrational energy levels have been examined. The results of these test calculations 
are listed in Table 2. First, we note that including all the potential energy terms and keeping 
the simple harmonic oscillator expression for the kinetic energy, cf. the column "Harm.+V", 
leads to an imbalance in the included contributions which results in an overestimation of 
the energy levels and in somewhat worse results than the simple harmonic approximation. 
Adding the anharmonic kinematic coefficients % spr and % sprt from Eq. (7) which, together 
with the potential energy function expanded up to and including the fifth order, lowers the 
energies making them closer to the experiment. Finally, it is interesting to examine the effect 
of adding the fifth-order terms Xsprtu in Eq. (5). It turns out that within the latter approxi- 
mation the errors are reduced only by 1-2 cm~' for the lowest energy levels. Therefore, 
we decided to skip Xsprtu m me production calculations. Including of the fifth-order kine- 
matic coefficients might be necessary for higher vibrational energies, especially in view of 
the fact that the absolute error for the series OV2O increases nonlinearly with Vo- However, 
when dealing with the expansion of the kinetic and potential terms, it is important to take 
all values of the same order in the Taylor expansion series, because neglecting even some of 
them may lead to high discrepancies in the final results. This follows from the fact that the 
variational procedure is very dependent on the non-diagonal terms included in the Hamil- 
tonian matrix. Finally, it is interesting to notice that Approximation 2a gives worse results 
than Approximation 2 for two high-/ cases presented in the Table, while for all other cases 
the latter approach gives slightly worse agreement with the experiment. 

For linear molecules the bending coordinates {Q2a^Qib} are present only as even pow- 
ers in the Taylor expansions, so V(Q2a, Qib) is even with respect to Qi a and Qib- This fact is 
due to the symmetry of the degenerate bending mode with respect to the plane perpendicular 
to the plane of vibration. We can also note an interesting property of high order force con- 
stants containing the bending coordinates, coming from the equivalence of the two bending 
vibration planes: 

-^2a,2 ti ,2„,2a ^2^,2^,2^,2^, ^2 a ,2 a ,2a,2fc/3. 

The pure vibrational energy levels relative to the ground-state energy level are listed in 
Tables 3-5 for the molecules HCN, OCS, and HCP, respectively. In these tables the funda- 
mental energy levels as well as some overtones are presented. Three approximations have 
been used in these tables: harmonic (Approximation 0), and Approximations 1 and 2 dif- 
fering in the way in which the V operator was truncated, cf. Section 6. The results are 
compared with the experimental data when available. The experimental data for the ener- 
gies and rotational constants were taken from Ref. [30] for HCN, Refs. [36, 79, 80] for OCS, 
and Refs. [52, 53] for HCR Additionally, Figures 1-3 show selected energy levels for lowest 
few bending modes of HCN, HCR and OCS calculated via all three methods and compared 
with the experiment. Differences between the theoretical and the experimental results are 
also reported for better readability of the tables. 

The results reported in Tables 3-5 show that in the series of approximations, starting 
with the harmonic and ending with Approximation 2, an increasing accuracy is observed in 
nearly all cases. Obviously, the harmonic approximation gives the largest errors with respect 
to the experimental values, and Approximation 1 always improves the agreement with the 
experiment. The level of improvement varies substantially depending on the molecule and 
the vibrational mode studied, but often the use of Approximation 1 causes even a five-fold 
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lowering of the energy gap between the calculated and experimental energy levels. Approx- 
imation 2 leads in almost all cases to results even closer to the experiment than Approxi- 
mation 1. It is also interesting to note that our calculations support the claim concerning the 
high anharmonicity of the C-H stretching mode. For the HCN and HCP molecules these 
modes are easily identified by looking at the Aq column. They are the only states with er- 
rors larger than 100 cm -1 . Although both Approximations 1 and especially 2 do a good 
job of bringing this mode close to the experiment, it could be envisaged that the next ap- 
proximation would give an even better agreement. Unfortunately, the analytical gradients 
within the CCSD(T) approach should be used for this purpose. Another way to improve 
the potential energy surface is to use Morse-like potential for this mode, which corresponds 
approximately to the C-H stretch. Such potential better describes upper vibrational states 
close to the dissociation limit, see for example Ref. [81]. Unfortunately, this leads to serious 
complications in the calculations of the matrix elements and has not been attempted in this 
paper. 

In Tables 3-5 a comparison of our results with previous calculations for every molecule 
is also presented. As already mentioned above, HCN is the most studied molecule and a vast 
theoretical material is available in the literature dealing with the problem of the vibrational 
HCN spectra. In Table 3 we listed values from four different sources. The first column with 
the theoretical results corresponds to the earliest work utilizing an ab initio PES [81], where 
one of the modifications of the DVR method was used to calculate the vibrational energies. 
The next set of results stems from a more recent implementation of the DVR method [82]. 
It should be stressed again that although the DVR method allows finding of vibrational 
energies up to very high excitation levels, it is limited to triatomic molecules and is therefore 
less general than the method presented in this paper in the curvilinear coordinates. Another 
work applying the DVR approach is Ref. [4], where DVR was implemented for the Watson 
Hamiltonian [3, 70, 76]. The last set of data was produced by Wang et al. [83] by applying 
a finite element method for the Watson Hamiltonian with J = 0. A perusal of this part of 
Table 3 allows us to conclude that the present vibrational energies are generally better than 
those of Ref. [81] and of the same level of accuracy than in the other references, except 
DVR implementation made in [82] where the accuracy is better, especially for combination 
bands, like (2 0° 3) in Table 3. This fact is reassuring in view of the just mentioned possible 
high-precision of the DVR method for the triatomic molecules. 

For carbonyl sulfide much fewer theoretical investigations are available. Since the OCS 
molecule is rigid, the model with the neglected bending mode can be used to study vibra- 
tionally excited stretching states. Such an approach has been pursued by Peterson et al. [45], 
The vibrational energies obtained in Ref. [45] are listed in Table 4. It is not surprising that 
the discrepancy between these results with the experimental values is larger than in the case 
of our approach, where the mixed bending- stretching terms have been taken into account. To 
our knowledge it was the only ab initio available potential energy surface for excited vibra- 
tional states. The next column of Table 4 contains the data of Ref. [84], where an algebraic 
model was used to analyze and interpret the experimental rovibrational spectra of small 
and medium-sized molecules. This method employs the Lie algebra techniques to obtain 
an effective Hamiltonian operator which describes rovibrational degrees of freedom. Since 
it does not use any ab initio results, but instead is solely based on the fitting of the spec- 
troscopic data, the results obtained there are very close to the experimental values. Finally, 
the last theoretical column in Table 4 contains data from Ref. [41], where a semi-classical 
model based on the adiabatic switching method (ASM) was used together with a spectro- 
scopic potential energy function. Also, in this case the good agreement with the experiment 
can be attributed to the use of some spectroscopic data in the actual calculations. Taking this 
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into account we can conclude that our method reproduces the vibrational energies for OCS 
much better that the other ab initio approach [45]. 

The DVR method was also applied to the phosphaethyne molecule [85]. The solution to 
the vibrational problem with an ab initio PES gave ,in general, a slightly better agreement 
with the experiment than our results. The following column in Table 4 shows the results of 
Ref. [56], where an ab initio PES was used to computed the rovibrational energy levels with 
a mixed variational-pertubation theory approach. Also here a somewhat better agreement 
between the vibrational energy levels and the experiment is found. However, the constants 
of the rovibrational interactions were calculated only for the ground state in this paper. 

Several possible reasons of the remaining discrepancies should be considered. The first 
issue is the quality of the CCSD(T)/cc-pVQZ electronic energies. Since the CCSD(T) method 
is known to be very accurate for molecules (unless very large distortions to the geometry are 
studied) and the cc-pVQZ orbital basis is also quite accurate, we can exclude the possibility 
that inaccurate electronic energies are the reason for the remaining errors. The most prob- 
able cause of the discrepancies are the truncation of the Taylor expansions which we have 
used for the g-tensor and for the potential V in the Hamiltonian. For instance, we used the 
second-order (i.e. linear and quadratic) expansion in the vibrational coordinates of [i a [s{Q), 
which may be insufficient for some excited rovibration bending energy levels, if the average 
geometry of the molecule differs considerably from the equilibrium geometry. 

For the C=N, C=P, and C=S stretching modes, generally larger errors are observed for 
the fundamental transition than for the other two. Errors of about 5 cm -1 can be seen in the 
former case and only 1-2 cm for the latter. The error increases approximately twofold for 
the first overtone (20°0). According to Koput et al. [56] this discrepancy can be explained by 
the effect of neglecting the correlation between the core electrons of the heavier atoms. This 
hypothesis is supported by the fact that only the Vi mode is so atrongly affected. It seems that 
further studies of the molecules containing atoms from the third row of the periodic table 
should account for the core correlation. Similar discrepancies have been found by Martin et 
al. [42]. 

Finally, one more complication arises from the fact that the values of the matrix elements 
p.™p are very sensitive to small non-diagonal anharmonic elements of the potential function 
V. These elements have a small influence on the calculated rovibrational energy levels, but 
cause a considerable change in the wave functions, and, hence in the average geometry of 
the molecule and its rotational energy levels if the molecule is vibrationally excited. 

In Tables 6-8 the constants for the rovibrational energies calculated within Approxi- 
mation 2 are presented. The experimental values for these molecules are also given, when 
available. These constants have been obtained by first computing rovibrational energy levels 
for various values of J and K, and then by performing fitting of the E r values to the ana- 
lytic form, cf. Eqs. (24)-(27). The J values up to and including 7=15 were used. The results 
show that the most important constants, i.e. B v , D v , and q v , are reproduced quite accurately. 
For the higher-order constants, the accuracy of the calculated rovibrational energies and the 
truncation of the maximum J was found to be insufficient, so the H v and q v j constants are 
not reported in Tables 6-8. Apart from the comparison to the experiment it is interesting 
to examine the dependence of the rotational constants on the vibrational quantum numbers. 
Such a dependence is especially visible for molecules containing a hydrogen atom, like HCN 
and HCP. The results allow us to conclude that the constants B v and D v grow linearly with 
V2, the quantum number for the degenerate bending mode, while for fixed Vo they decrease 
with |/|, the vibrational angular momentum quantum number for the bending mode Vt - This 
finding can be explained by the fact that for a molecule in an excited bending mode, the 
component of I : becomes smaller which leads to a lowering of the constants B v and D v . On 
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the other hand, decreasing the constants while increasing |/| is mostly caused by the interac- 
tion of the total angular momentum with the vibrational one. These constants also become 
smaller when the quantum numbers Vi and V3 increase. This result can be rationalized by 
the fact as the vibrationally averaged interatomic distance (r) becomes larger, the distance 
differs more from the equilibrium distance. This assumption can be supported by calculating 
the vibrationally averaged values of the bond distances fo the molecules under study. Luck- 
ily, such calculations were reported by Laurie et al. [86] who presented (rue) = 1 -0739, 
(>"cn) = L1574 A for the ground state. These values are larger than the equilibrium values 
from Table 1 . 



8 Summary 

The rovibrational Hamiltonian in curvilinear vibrational coordinates has been used to solve 
the nuclear motion problem for several linear triatomic molecules in curvilinear vibrational 
coordinates. The curvilinear coordinates were shown to have advantages over the more com- 
monly used linear coordinates. Comparison with the exising theoretical results shows that 
the present approach works equally well comparing to other ab initio calculations of vi- 
brational energy levels, and is not limited to triatomic molecules. In its current form the 
approach presented in this paper provides a good accuracy for most vibrational levels and 
recovers the rotational and the /-doubling constants for the studied molecules. Further im- 
provements to the code will also include the calculation of intensities, so that the full spec- 
trum of the molecules can be simulated. 
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Table 1 Equilibrium bond lengths (in A) obtained in this work and taken from the experiment. r\ (ri) denotes 
the distance between carbon and a lighter (heavier) atom. 



Molecule 


''l 


T2 


Reference 


HCN 


1.06686 


1.15646 


this work 




1.065 


1.153 


[27] 


HCP 


1.07216 


1.54732 


this work 




1.070 


1.540 


[54] 


OCS 


1.15830 


1.56902 


this work 




1.156 


1.561 


[44] 



Table 2 The influence of the kinematic coefficients on the vibrational energies for the HCN molecule. All 
values in cm~ 1 . Approximation 2a with added Xsprtu 



Vl VjV3 


Harm. 


Harm. + V 


Appr. 2 


Appr. 2a 


Exp. 


1' 


721.9 


727.6 


713.4 


712.9 


712.0 


2°0 


1443.9 


1447.9 


1415.4 


1413.7 


1411.4 


2 2 


1443.9 


1460.4 


1424.5 


1423.3 


1426.1 


3 1 


2165.8 


2170.4 


2119.6 


2117.7 


2113.5 


3 3 


2165.8 


2190.8 


2145.6 


2143.3 


2157.2 


4°0 


2887.8 


2926.5 


2820.1 


2816.1 


2803.0 


0° 1 


2123.0 


2148.5 


2098.1 


2098.0 


2096.4 


1 0°0 


3435.3 


3554.5 


3316.8 


3314.6 


3311.4 


0° 2 


4246.0 


4231.5 


4168.2 


4165.3 


4173.2 


2 0°0 


6870.7 


6934.7 


6529.4 


6524.9 


6519.5 


l 1 1 


2844.9 


2898.3 


2810.5 


2807.4 


2805.6 


2 0° 3 


13239.6 


13432.9 


12675.1 


12669.2 


12658.0 



21) 



Table 3 Experimental and theoretical vibrational energy levels for the HCN molecule. Theoretical results 
were obtained within the harmonic approximation (Approximation 0), and with Approximation 1 and Ap- 
proximation 2 defined in the text. A; denotes the error of a given approximation ; with respect to the exper- 
iment, i.e. Aj = £(Appr. i) — £(Exp). All values are in cm -1 . Experimental data are taken from Ref. Ref. 
[30]. 



Vi^v? Harm. A Appr. 1 A { Appr. 2 A 2 Exp Ref. [ 8 1 1 Ref. [82] Rel. [4| Ki'l.|~ 



1' 


721.9 


9.9 


715.2 


3.2 


713.4 


1.4 


712.0 


718.4 


712.0 


715.9 


713.0 


2°0 


1443.9 


32.5 


1421.1 


9.7 


1415.4 


4.0 


1411.4 


1418.9 


1411.4 


1414.9 


1421.6 


2 2 


1443.9 


17. S 


1422.1 


-4 


1424.5 


-1.6 


1426.1 


n/a 


n/a 


n/a 


n/a 


3 1 


2165.8 


52.3 


2125.6 


12.1 


2119.6 


6.1 


2113.5 


2126.0 


2113.45 


n/a 


2127.9 


3 3 


2165.8 


21.8 


2135.5 


-8.5 


2145.6 


1.6 


2157.2 


n/a 


n/a 


n/a 


n/a 


4°0 


2887.8 


85.8 


2850.4 


37.4 


2820.1 


18.1 


2803.0 


2812.5 


2807.1 


2801.5 


2834.3 


0° 1 


2123.0 


26.6 


2095.1 


-1.3 


2098.1 


1.7 


2096.4 


2090.3 


2091.0 


2100.6 


2083.2 


I0°0 


3435.3 


123.9 


3313.4 


2 


3316.8 


5.4 


3311.4 


3334.1 


3311.5 


3307.7 


3340.4 


0° 2 


4246.0 


72.8 


4181.7 


8.5 


4168.2 


-5 


4173.2 


4161.5 


4173.1 


4181.5 


4146.3 


20°0 


6870.7 


351.2 


6598.4 


78.9 


6529.4 


-9.9 


6519.5 


6553.2 


6519.6 


6513.5 


a/a 


1 1 1 


2844.9 


39.3 


2809.4 


3.8 


2810.5 


4.9 


2805.6 


2806.4 


2807.1 


n/a 


a/a 


2 0° 3 


13239.6 


581.6 


12795.1 


137.1 


12675.1 


17.1 


12658.0 


n/a 


12658.0 


n/a 


a/a 


2° 1 


3566.9 


64.8 


3532.1 


30.0 


3512.0 


9.9 


3502.1 


3507.4 


3511.0 


3511.0 


3511.0 


2 2 1 


3566.9 


44.2 


3541.7 


19.0 


3529.4 


6.7 


3522.7 


n/a 


n/;i 


n/a 


a/a 


1 3' 


5601.1 


232.8 


5435.7 


67.4 


5398.4 


30. 1 


5368.3 


5387.1 


n/a 


n/a 


a/a 



Table 4 Experimental and theoretical vibrational energy levels for the OCS molecule. All values are in cm 
See Table 3 for an explanation of the symbols. Experimental data are taken from Ref. [79]. 



V1V2V3 


Harm. 


4, 


Appr. 1 


42 


Appr. 2 


4 3 


Exp 


Ref. [45] 


Ref. [84] 


Ref. [41] 


1 1 


524.2 


3.8 


522.1 


1.7 


520.9 


0.5 


520.4 


n/a 


n/a 


520.3 


2°0 


1048.4 


1.4 


1047.3 


(1.3 


1047.1 


0.1 


1047.0 


n/a 


1047.0 


1046.9 


2 2 


1048.4 




1045.7 




1042.1 




n/a 


n/a 


n/a 


1048.8 


3' 


1572.7 


-0.7 


1573.1 


-0.3 


1573.7 


0.3 


1573.2 


n/a 


1569.0 


1573.4 


3 3 


1572.7 




1573.4 




1574.0 




n/a 


n/a 


n/a 


1561.5 


04° 


2096.9 


-7.9 


2101.3 


-3.5 


2102.4 


-2.4 


2104.8 


n/a 


2088.0 


2105.8 


0° 1 


871.7 


12.4 


861.4 


2.1 


859.1 


-0.2 


859.3 


869.0 


855.5 


859.2 


1 0°0 


2094.7 


32.4 


2078.4 


16.2 


2067.1 


4.9 


2062.2 


2071 


2057.9 


2062.1 


0° 2 


1743.4 


32.7 


1715.7 


5 


1712.7 


2 


1710.7 


1730.0 


1705.0 


1710.6 


2 0°0 


4189.3 


86.9 


4127.2 


25.8 


4110.7 


9.3 


4101.4 


4119.0 


4097.1 


4101.4 


1 1 1 


1395.9 




1379.4 




1373.1 


n/a 


n/a 


n/a 


n/a 


1372.8 


2 0° 3 


6804.2 




6704.1 




6687.5 




n/a 


n/a 


n/a 


n/a 


2° 1 


1920.1 


27.0 


1915.2 


23.0 


1899.7 


7.5 


1892.2 


n/a 


1895.5 


1891.7 


2 2 1 


1920.1 




1911.4 




1890.5 




n/a 


n/a 


1886.7 


1890.5 


1 3 1 


3667.4 


52.0 


3643.2 


27.8 


3618.8 


3.4 


3615.4 


n/a 


3617.9 


3614.4 



Table 5 Experimental and theoretical vibrational energy levels for the HCP molecule. All values are in cm 
See Table 3 for an explanation of the symbols. Experimental data are taken from Refs. [52]" and [53]*. 



V1V2V3 


Harm. 


Ao 


Appr. 1 


A, 


Appr. 2 


A 2 


Exp 


Ref. [85] 


Ref. [56] 


Ref. [84] 


1' « 


686.9 


11.95 


675.2 


0.2 


674.4 


-0.6 


675 


675 


675.2 


n/a 


2°0 " 


1373.9 


37.9 


1349.1 


13.1 


1340.4 


4.4 


1 336 


1 336 


1335.0 


1332.9 


2 2 " 


1373.9 


9.9 


1369.1 


5.1 


1365.7 


1.7 


1 364 


1364 


n/a 


1340.2 


3'0 « 


2060.8 


58.8 


2020.4 


18.4 


2008.7 


6.7 


2002 


2002 


2002.1 


1995.6 


3 3 " 


2060.8 




2041.1 




2036.1 




n/a 


n/a 


n/a 


n/a 


4° " 


2747.8 


105.0 


2721.1 


78.3 


2715.8 


73.0 


2643 


n/a 


2653.8 


n/a 


0° 1 " 


1294.9 


15.9 


1287.2 


8.2 


1285.1 


6.1 


1279 


n/a 


1275.5 


1279.0 


1 0° " 


3345.3 


129.3 


3301.7 


85.7 


3207.1 


-8.9 


3217 


n/a 


3215.3 


n/a 


0° 2 " 


2589.8 


40.S 


2571.4 


22.4 


2570.1 


21.1 


2549 


n/a 


2540.1 


2547.4 


2 0° « 


6690.5 




6671.4 




6667.8 




n/a 


n/a 


6321.0 


n/a 


1 1 1 " 


1981.8 


34.S 


1960.1 


13.1 


1955.1 


8.1 


1947 


n/a 


n/a 


n/a 


2 0° 3 « 


10575.2 




10384.7 




10276.7 




n/a 


n/a 


n/a 


n/a 


2° 1 " 


2668.8 


70.S 


2634.2 


36.2 


2608.4 


10.4 


2598 


2598 


n/a 


2597.2 


2 2 1 " 


2668.8 


42.8 


2647.2 


21.2 


2636.3 


10.3 


2626 


2626 


n/a 


2604.5 


1 3 1 « 


5406. 1 




5302.7 




5249.7 




n/a 


n/a 


n/a 


n/a 
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Table 6 Rovibrational constants for HCN. Theoretical results were obtained with Approximation 2. All val- 
ues are in cm . Experimental data are taken from Ref. [30]. 





B v Exp 


By Calc 


D v x 10" Exp 


D v x 10" Calc 


q v x 10 3 Exp 


q Y x 10 J Calc 


0" 


1.478222 


1.479782 


2.91028 


2.59259 






l 1 


1.481773 


1.481931 


2.97746 


3.01327 


7.48773 


7.51431 


2°0 


1.485828 


1.487207 


3.04799 


3.08957 


7.59561 


7.63292 


2 2 


1.484997 


1.483982 


3.04079 


3.11047 






03 1 


1.489575 


1.492028 


3.11365 


3.19783 


7.70926 


7.81492 


3 3 


1.487868 


1.491873 


3.10024 


3.29147 






l 1 1 


1.471574 


1.475982 


2.98254 


3.11873 


7.47954 


7.49425 


2° 1 


1.475493 


1.482542 


3.05017 


3.78874 


7.57341 


7.61422 


2 2 1 


1.474678 


1.491012 


3.04520 


3.85578 






1 3 1 


1.479809 


1.481287 


3.11952 


4.50034 


7.85215 


7.83421 



Table 7 Rovibrational constants for OCS. Theoretical results were obtained with Approximation 2. All val- 
ues are in cm~'. Experimental data are taken from Refs. [80]° and [36]*. 



VlViV 3 




By Exp 


B v Calc 


D v x 10 s Exp 


D v x 10 s Calc 


q v x 10 4 Exp q v : 


< 10 4 Calc 


0" 


a 


0.202856 


0.202506 


4.341064 


4.703412 






l 1 


a 


0.203209 


0.203347 


4.411467 


5.198434 


2.12193 


2.11932 


2°0 


b 


0.203480 


0.203135 


3.66256 


3.963137 


n/a 


2.13421 


2 2 


b 


0.203559 


0.203397 


5.23854 


4.598443 






3 1 




n/a 


0.203492 


n/a 


4.993237 


n/a 


2.16402 


3 3 




n/a 


0.204431 


n/a 


5.194394 






l 1 1 


a 


0.202657 


0.203054 


4.54206 


5.94217 


2.28496 


2.25421 


2° 1 


a 


0.202953 


0.203384 


4.55522 


5.53215 


2.2238 


2.18401 


2 2 1 


a 


0.203048 


0.203043 


4.63515 


5.73263 






1 3 1 




n/a 


0.203349 


n/a 


6.143694 


n/a 


2.17754 



Table 8 Rovibrational constants for HCR Theoretical results were obtained with Approximation 2. All values 
are in cm -1 . Experimental data are taken from Ref. [53]. 





By Exp 


B r Calc 


Dy x 10" Exp 


Dy x 10" Calc 


q r x 10 3 Exp 


q r x 10 J Calc 


0" 


0.666326 


0.667483 


0.70392 


0.93820 






l 1 


0.666771 


0.667333 


0.71196 


0.99482 


1.63032 


1.65339 


2°0 


0.667344 


0.668391 


0.719160 


1.03945 


1.65970 


1.68432 


2 2 


n/a 


0.663641 


n/a 


1.09384 






3 1 


0.667784 


0.668380 


0.72502 


1.10939 


1.70887 


1.73283 


3 3 


n/a 


0.673452 


n/a 


1.03988 






l 1 1 


0.662932 


0.665482 


0.72052 


0.928438 


1.56214 


1.49844 


2° 1 


n/a 


0.668391 


n/a 


0.91938 


n/a 


1.61925 


2 2 1 


n/a 


0.663942 


n/a 


1.39492 






1 3 1 


n/a 


0.653393 


n/a 


1.03821 


n/a 


1.66754 
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